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Abstract 

The problem of Laplacian growth in two dimensions is considered within the Loewner-equation 
framework. Initially the problem of fingered growth recently discussed by Gubiec and Szymczak 
[T. Gubiec and P. Szymczak, Phys. Rev. E 77, 041602 (2008)] is revisited and a new exact solution 
for a three-finger configuration is reported. Then a general class of growth models for an interface 
growing in the upper-half plane is introduced and the corresponding Loewner equation for the 
problem is derived. Several examples are given including interfaces with one or more tips as well as 
multiple growing interfaces. A generalization of our interface growth model in terms of "Loewner 
domains," where the growth rule is specified by a time evolving measure, is briefly discussed. 
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I. INTRODUCTION 

The Loewner equation [1] is an important result in the theory of univalent functions 
that has found important applications in nonlinear dynamics, statistical physics, and 
conformal field theory {^-[(j]. In its most basic formulation, the Loewner equation is a first- 
order differential equation for the conformal mapping gt{z) from a given "physical domain," 
consisting of a complex region P minus a curve T t emanating from its boundary, onto a 
"mathematical domain" represented by P itself. Usually, P is either the upper half-plane or 
the exterior of the unit circle, but recently the Loewner equation for the channel geometry 
was also considered 7J. The Loewner equation depends on a driving function, here called 
a(t), that is the image of the growing tip under the mapping gt(z). An important develop- 
ment on the theory of the Loewner equation was the discovery by Schramm that when 
the driving function a(t) is a Brownian motion the resulting Loewner evolution describes 
the scaling limit of certain statistical mechanics models. This result spurred great interest 

nn 

in the so-called stochastic Loewner equation [4-6]. 

Recently, the deterministic Loewner equation was used to study the problem of Laplacian 



fingered growth in both the half-plane and radial geometries P.lloj| as we 11 as in the channel 
geometry 7|. In this class of models the growth takes place only at the tips of slit-like 
fingers and the driving function a(t) has to follow a specific time evolution in order to 
ensure that the tip grows along gradient lines of the corresponding Laplacian field. In spite 
of its simplicity, the model was able to reproduce some of the qualitative behavior seen in 
experiments on fingered growth 0, • Nevertheless, treating the fingers as infinitesimally 
thin is a rather severe approximation and so this "thin" finger model is applicable only if the 
width of the fingers is much smaller than the separation between fingers. On the other hand, 
there are other instances of Laplacian growth where "extended fingers" (i.e., fingers with a 
finite, non-negligible width) are observed [12|] and for such description in terms of 

Loewner evolutions is still lacking. 

One of the main motivations of this paper is to seek to develop a framework based on the 
Loewner equation in which one can study more general growth problems in two dimensions. 
We start our analysis by first revisiting the problem of fingered growth discussed in [?| and 
report here a novel exact solution for the case of a symmetrical configuration with three 
fingers. We then move on to discuss a general class of growth models where an interface 
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grows into the upper-half plane starting from a segment on the real axis. In our model, the 
growth rate is specified at a certain number of special points along the interface, referred to 
as "tips" and "troughs", which then determine the growth rate at the other points of the 
interface according to a specific growth rule (formulated in terms of a polygonal curve in the 
mathematical plane). By making appropriate use of the Schwarz-Christoffel transformation, 
we are then able to derive the Loewner equation for the problem. Several examples are 
given in which the interface evolution is numerically computed by a direct integration of 
the Loewner equation. We also briefly discuss a more general formulation of the Loewner 
equation for interface growth in the upper-half plane where the growth rule is given in terms 
of a time evolving measure. This approach may open up the possibility for studying other 
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interesting growth problems, such as the growth of fractal interfaces 
of stochastic Loewner evolutions. 

The paper is organized as follows. In Sec. [II] we briefly review the problem of fingered 
growth in the upper half-plane in the context of the chordal Loewner equation and report 
a new solution for a symmetrical configuration with three fingers. In Sec. II I II we discuss a 
general class of growth models in which the growing domain is delimited by an interface 
in the upper half-plane. Here we first derive in Sec. IIII Al a new Loewner equation for the 
case in which the growth rule is specified in terms of a polygonal curve in the mathematical 
plane. Several several examples are given in Sec. IIII B\ including the cases of interfaces with 
one or more tips as well as the case of multiple growing interfaces. We then briefly discuss in 
Sec. IIII CI a more general formulation for interface growth in the upper-half plane in terms 
of "Loewner domains." Our main results and conclusions are summarized in Sec. IIVI 



II. FINGERED GROWTH AND THE CHORDAL LOEWNER EQUATION 

In order to set the stage for the remainder of the paper and to establish the relevant 
notation, we begin our discussion by briefly reviewing the problem of fingered growth in the 
context of the chordal Loewner equation. To this end we consider first the simplest Loewner 
evolution, namely, that in which a curve starts from the real axis at t — and then grows 
into the upper half-z-plane HI, where 

H = {z = x + iy e C : y > 0}. 
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FIG. 1: The physical z-plane and the mathematical w- and ^-planes at times t and i+r, respectively, 
for a single finger in the upper half-plane. The mapping gt maps the curve Ft onto a segment of 
the real axis on the w-plane, whereas the accrued portion of the curve during the infinitesimal time 
interval r is mapped to a vertical slit. The mapping gt is obtained as the composition of gt+ T an d 
the slit mapping F; see text. 

The curve at time t is denoted by T t and its growing tip is labeled by j(t). Now let gt(z) 
be the conformal mapping that maps the "physical domain," corresponding to the upper 
half-z-plane minus the curve T t , onto the upper half-plane of an auxiliary complex w-plane, 
called the "mathematical plane," i.e., we have w = gt(z), where 

g t : M\T t -> H, (1) 

with the curve tip 7(t) being mapped to a point a(t) on the real axis in the w-plane, as 
shown in Fig. [TJ Furthermore, we consider the growth process to be such that the accrued 
portion of the curve from t to t + r, where r is an infinitesimal time interval, is mapped 
under gt{z) to a vertical slit in the mathematical w-plane; see Fig. [TJ The mapping function 
gt(z) must also satisfy the initial condition 

go(z) = z, (2) 

since we start with an empty upper half-plane. We also impose the so-called hydrodynamic 
normalization condition at infinity: 

g t (z) = z + 0( T -r), z^oo. (3) 

These conditions specify uniquely the mapping function gt(z). 

From a more physical viewpoint, the problem formulated above belong to the class of 
Laplacian growth models where an interface evolves between two phases driven by a scalar 
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field (f>(x, y; t), representing, for example, temperature, pressure, or concentration, depending 
on the physical problem at hand. In one phase, initially occupying the entire upper half- 
plane, the scalar field <fi satisfies the Laplace equation 

V 2 = 0, (4) 

whereas in the other phase one considers </>=const., say = 0, with the curve T t representing 
a finger-like advancing interface between the two phases. (Here the finger is assumed to 
be infinitesimally thin.) The complex potential for the problem can then be defined as 
w(z,t) = i/}(x,y;t) + i<p(x,y;t), where ip is the function harmonically conjugated to 4>. On 
the boundary of the physical domain, consisting here of the real axis together with the curve 
r t , we impose the condition = 0, whereas at infinity we assume a uniform gradient field, 
V0 ~ y, or alternatively, 

w(z,t) fa z, z — > oo. (5) 

From this point of view, the mapping function gt(z) introduced above corresponds precisely 
to the complex potential w(z,t) of the problem. In particular, the fact that in the w-plane 
the curve grows along a vertical line implies that the finger tip grows along gradient lines 
in the z-plane. To specify completely a given physical model, one has also to prescribe the 
interface velocity, which is usually taken to be proportional to some power rj of the gradient 
field: 

Vn ~ (6) 

We anticipate that for a single growing finger the specific velocity model is not relevant in 
the sense that the finger shape will be independent of the exponent r], which only affects 
the time scale of the problem. (For multifingers, however, different 77's may yield different 
patterns 3].) 

For convenience of notation, we shall represent the mathematical plane at time t + r as 
the complex ^-plane and so we write £ = g t + T (z). Now consider the mapping w = F(Q, 
from the upper half-^-plane onto the mathematical domain in the w-plane; see Fig. [TJ The 
mapping function g t+T (z) can then be given in terms of gt(z) by 

g HT = F- 1 og t , (7) 

where F" 1 is the inverse of F((). The above relation governs the time evolution of the 
function gt(z) and naturally leads to the Loewner equation. A standard way of showing this 
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is to construct the slit mapping F(Q explicitly, substitute its inverse in ([7]), and then take 
the limit r — > 0. One then finds the so-called chordal Loewner equation: 

9t{z) = dit) (8) 
g t (z) - a(t) 

together with the condition 

a(t) = 0, (9) 

so that a(t) = a = const., which implies that the tip jit) = g^i^ait)) simply traces out a 
vertical line in the z-plane. The growth factor d(t) is related to the tip velocity by [7| 

d(t) = \ti'(a(t))r /2 -\ (10) 

where ft(w) is the inverse of gt(z). As already anticipated, we can take d(t) = d Q (which 
corresponds to rj = —2) by reparametrizingthe time coordinate. 

In the case of multiple growing fingers |7j, the chordal Loewner equation becomes 

* = £-^> (H) 
with the time evolution of the singularities cii(t) given by 

= (12) 

We recall that conditions ( JT2l) ensure that the fingers grow along gradient lines, meaning 
that the path traced by each finger tip 7i(t) from time t to t + r is mapped by gt{z) onto a 
vertical slit. Notice also that if the growth factors di(t) are all the same they can be assumed 
to be constant, say, di(t) = 1 (which corresponds to 77 = —2), since this amounts to a mere 
rescaling of the time variable. For two symmetrical fingers, i.e., di = <i 2 , equation ffTTj) can 
be integrated exactly to yield the mapping function gt(z), from which the finger shapes can 
be computed analytically [2]. A related exact solution for two fingers was reported in [jjj]. 

An exact solution for gt{z) can also be obtained for a symmetrical configuration with three 
fingers where the middle finger grows vertically, say, along the y axis, and the two flanking 
fingers are the mirror images of one another with respect to the y axis. This situation is 
attained by choosing identical growth factors, ^1 = ^2 = ^3 = ^0; and symmetrical initial 
conditions, 03(0) = — ai(0) = qq and 02(0) = 0. This initial symmetry is, of course, preserved 
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by the dynamics and so we have a,2{t) = and a^(t) = —ai(t) = a(t) for all t > 0, which 
implies from ffT2]) that 



a(t) = ^Jal + 3d t. (13) 

Using this result one can integrate exactly the Loewner equation ( ITT]) and obtain, after a 
straightforward calculation, the positions 7^3 (t) of the tips of the flanking fingers in terms 
of the following implicit equation 

hiM - P + \ 1+a biM - ^Y~ a = (i - - ^f- a {al + 3d tf, (i4) 

where a = ^3/19 and ± = (9 ± v / 57)/6. In flJJD one considers only the roots that lie 



in the upper-half plane and satisfy the initial condition 713 = =Fao. In Fig. 2(a) we plot 
this solution for ao = 1 and do — 1. From this figure one can see that the side fingers are 
'repelled' by the middle finger and tend to straight lines for sufficiently large times. In fact, 
one can verify from (}T4"j) that in the asymptotic limit t — > 00 the flanking fingers approach 
the straight lines arg(z) = f (1 + y/3/19) and arg(z) = |(3 - y/3/19). 

If the growth factors are kept all equal but the initial condition is no longer symmetric, 
then the middle finger will initially be more repelled by the closest finger, thus leading to 
a symmetrical configuration in the limit t — > 00 where the asymptotic (vertical) position of 
the middle finger now depends on the initial condition. On the other hand, if the growth 
factors are not all the same then the finger with greatest di will grow faster and screen 
the others. Here, however, the screening is partial in the sense that the ratio between the 
velocity of a slower finger and that of the fastest finger reaches a positive constant (different 



from unity). An example of this case is shown in Fig. 2(b) for the early stages of a three- 
finger evolution. In order to generate the curves shown in this figure we used the numerical 
scheme described in [yj]. More specifically, we start with the "terminal condition" g t = a>i(t) 
and integrate the Loewner equation ffTTl) backwards in time, using a Runge-Kutta method 
of second order, to get the position of the respective tip 7,;(f) = g . We dealt with the 
pole singularities in (ITT]) in a manner similar to that used in [14j, which in our case meant 
finding the appropriate root of a fourth-degree polynomial that gives the value of g at time 
t — dt, from where the Loewner can be integrated backwards to yield go- It is interesting to 



notice that the curves in Fig. 2(b) resemble certain fingering patterns seen in combustion 



experiments, see, e.g., Fig. le of Ref. jllj . even though the solutions above are valid for the 



upper-half plane, whereas the experiments take place in a channel geometry, and we have 
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FIG. 2: Three fingers growing in the upper-half plane. In (a) we have a symmetrical configuration 
with d\ = d,2 = d>2 = 1 and initial positions of the fingers at 03(0) = — ai(0) = 1 and 02(0) = 0; 
in (b) we have d\ = 1, &i = 5, d% = 0.9, while the initial positions are ^3(0) = —a\(0) = 1 and 
a 2 (0) = -0.4. 

used constant growth factors (although with different values for each finger) rather than 
equation f lTOj) . 

Before leaving this section we wish to recall here that the evolution entailed by the chordal 
Loewner equation, as given by flB]) or f JTTj) . describes a rather general class of "local growth" 
processes [5[ in which there is one curve or a set of curves growing from the real axis into 
the upper-half plane EI. Indeed, it is a theorem js, 6| that if T^oo] is a simple curve from 
the origin to 00 in the upper-half plane H, then there exists a function a(t) such that the 
curve is generated by a Loewner evolution. Conversely, if a(t) is a sufficiently well-behaved 
function (more specifically, Holder continuous with exponent > 1/2), then j t = g^i^ait)) is 
a growing curve in HI. More general growth processes not restricted to growing curves can 
be generated by the so-called "Loewner chains" jsj: 

p(x)dx 



<nr.) -- I —r- , (15) 

g t (z) - x 

where the density of singularities p(x) can be viewed as a measure of the growth rate at a 
point z at the boundary of the growing set that is the preimage of x under gt(z). Since the 
shape of the growing domain is fully encoded in the map gt(z), equation ( jl~5i) specifies the 
growth model once the density p(x) is known. Although the formalism of Loewner chains 
is rather general, its practical usefulness is somewhat limited by the requirement that the 
relevant density p(x) must be specified a priori. In this context, it should be noted that the 
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only specific Loewner evolution that has been studied in detail concerns the local growth 
models mentioned above where the density is a sum of Dirac 5 functions ^J. In the next 
section, we will study a new class of interface growth models where the growth velocity is 
specified at certain points of the interface, which then determine the growth rate at all the 
other points according to a specific rule. 



III. INTERFACE GROWTH IN THE HALF-PLANE 

A. Generalized Loewner equation for a growing interface 

Here we consider the problem of an interface starting initially from a segment, say, [—1, 1], 
along the real axis and growing into the upper half-z-plane, as indicated in Fig. [3j We 
suppose that the growing interface has a certain number of special points, referred to as 
tips or troughs, where the growth rate is a local maximum or a local minimum, respectively, 
while the interface endpoints remain fixed, as illustrated in Fig. [3] for the case of an interface 
with two tips (points B and D) and one trough (point C). Let us now denote by T t the 
interface at time t and by K t the growing region delimited by T t and the real axis. Here we 
assume that the curve T, is simple so that the domaiu D, = U\K, is simply counected. (In 
more technical terms, K t is a hull [5|.) We then consider the mapping w = gt(z) from the 
physical domain D t in the z-plane to the upper-half plane EI in the mathematical u>-plane: 

where the mapping function gt{z) is required to satisfy the hydrodynamic normalization 
condition 

g t (z) = z + 0(^-), z^oo, (16) 
together with the initial condition 

9o(z) = z. (17) 

By definition, the interface T t is mapped under gt{z) to an interval [cii(t), aiv(t)] on the real 
axis of the w-plane, where ai(t) and ajv(i) are the images of the endpoints z = ±1, while 
the tips and troughs are mapped to the points a,j(t), i = 2, N — 1, with N being the total 
number of special points (i.e., tips, troughs and endpoints) of the interface. For instance, in 
Fig. [3] we have N = 5. 
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FIG. 3: The physical and mathematical planes for a growing interface. The mapping gt(z) maps 
the interface Tt to a segment on the real axis, while the interface r\+ T is mapped to a polygonal 
curve; see text for details. 

The growth dynamics is specified by requiring that the tips and troughs of the interface 
grow along gradient lines, while the interface endpoints remain "pinned" at z = ±1, in such 
a way that the interface T t+T at time t + r, for infinitesimal r, is mapped under gt{z) to a 
polygonal curve in the w-plane, as shown in Fig. [31 The domain D t+r = M\K t+T is mapped 
under gt(z) to a degenerate polygon whose interior angle at the z-th vertex is denoted by 
7r(l — ctj), with the convention that if the angle is greater than 7r the corresponding parameter 
ati is negative. It is easy to verify that the parameters cVs satisfy the following relation 

N 

X)«< = °: ( 18 ) 

i=l 

where we recall iV is the number of vertices of the polygonal curve in the w-plane. If we 
now denote by hi the height of the z-th vertex as measured from the real axis (see Fig. [3]), 
it is clear that hi and hence the angle parameters a, all go to zero as r — > 0. In this limit 
one obtains another relation among the cVj's: 

N 

a,aj = 0. (19) 



i=i 



Let us know derive the Loewner evolution for the growth process described above. Since 
the domain in the w-plane has a polygonal s hap e, the mapping w = F(Q can be obtained 
from the Schwarz-Christoffel transformation 16]: 



9t = F(g. 



t+Tj 



gt+r 
Co 



2V 



Hi(- ai (t + T)]- a >d( + F(( ), 



(20) 



i=l 
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where Co is a point to be chosen appropriately; see below. The integral in ( 120]) cannot be 
performed exactly for arbitrary a^s, and so in order to obtain the Loewner equation for 
this case we need an alternative approach. First we expand the integrand in (120 p up to first 
order in the infinitesimal parameters ctj's, thus getting 

9t « J' t+r 1 1 - a < ln K - + r )} | d ( + F (Co), (21) 
which after integration becomes 

9t = 9t + r - Co - G(g t+T ) + G(Co) + F(Co), (22) 

where 

N 

G(C) =X;«*[C-Oi(« + r)] ln[C - a,(t + r)] (23) 

i=l 

Now expanding (1221 up to first order in r, taking r — >• 0, and using the boundary condition 
lim^^oo^ = 0, which follows from ([3]), one then obtains the following generalized Loewner 
equation 

N 

g t (z) = J2 d ^)i9t ~ Oi(t)) ln[g t - ai (t)}, (24) 

i=l 

together with the generic condition 

, im C.-F(Co)-G(C.) = (25) 
r^o r 

where the growth factors di(t) appearing in (1241 are defined by 

ry ■ 

dAt) = lim — . (26) 

r->0 T 

The equations governing the time evolution for the functions flj(t) can now be obtained 
by considering Co = (k{t + t) in ( J25l) . for i = 1, N, with F(Co) then being the respective 
image of a«(t + r) in the w-plane: 

F(a i (t + r)) = a i (t)+ih u (27) 

where for the endpoints we have h\ = = 0. Using ( 1271) in ( )25|) and performing a 
straightforward calculation, one finds 

hi = dj(t)(a,i — dj) In | a, — aj|. (28) 

3=1 
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In view of (126]) . equations ( ITS]) and (1X9]) for the parameters imply equivalent relations 
for the growth factors di(t): 

N 

£> = > (29) 

i=l 
N 

Mi = 0. (30) 

i=i 

From these equations we can then express the parameter functions, di(t) and d^(t), of the 
endpoints in terms of the other functions di(t), so that the growth model described by the 
Loewner evolution (124"]) and (|28|) is completely specified by prescribing the growth factors 
di(t) for the tips and troughs of the interface. [We remark parenthetically that although our 
growth model naturally allows for vertices that correspond to neither tips nor troughs we 
shall not consider such cases here; see however Sec. IIII O for a discussion about more general 
growth rules.] The parameter functions di(t) could in principle be related to the velocities 
of the corresponding points on the interface by equations analogous to ffTOl) . For instance, in 
the simplest case 77 = —2 we have di(t) = do- For general rj, however, the function di(t) has a 
complicated implicit dependence on the function g t (z) involving the second derivative of its 
inverse, which renders the numerical integration of the Loewner equation (1241) very difficult 
(if possible at all). On the same token, the alternative approach for computing Loewner 
evolutions based on direct iteration of slit mappings 0, [nj is unlikely to result very useful 
in the case of growing interfaces, because the relevant mapping to be iterated is not known 
in closed form; see (1201) . For these reasons, we shall consider in the examples below only 
the case di(t) = di = const, (although not all d^s need be the same), for which the Loewner 
equation ( 124"1) can be easily integrated. 



B. Examples 

We start by considering the case in which the interface has only one tip so that the 
number of vertices is obviously N = 3. We can then solve (129]) and (130]) to obtain the 
growth factors di(t) and d 3 (t) as a function of <i 2 (t): 

di(t) = _ Q 3^)-«2(t) = _ a 2 (t)-a 1 (t) 

U ' a 3 (t)- ai (t) 21 ; ' 31 ; a 3 (t)- ai {t) 2K) K J 

As mentioned above, the growth factor d,2(t) can in principle be related to the tip velocity 

but here the specific form of ^(t) is not relevant for the shape evolution, for it only defines 
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FIG. 4: Loewner evolution for a non-symmetric interface growing in the upper half-plane. The 
solid curves show the interface at various times t, starting from t = 0.5 up to t = 2.0, with a 
time separation of At = 0.3 between successive curves, whereas the dashed line indicates the tip 
trajectory. 

the time parameter, and so we set d,2(t) = — 1- (Recall that according to our convention the 
growth factors for tips are negative while those for troughs are positive.) As for the initial 
conditions, we have ai(0) = —1 and 03(0) = 1, corresponding to endpoints z = ±1, so that 
we have to specify only the initial location, 02(0), of the tip. If we start with the tip at the 
center, i.e., 02(0) = 0, we obtain a symmetrically growing interface [3], since in this case 
equations ( 1281) imply that 02 (t) = and a\(t) = —a^t) for all t. If, on the other hand, we 
start with the tip off-center, i.e., 02(0) 7^ 0, we then get an asymmetric growing interface 
(which can be regarded as an "extended finger" in the sense that it encloses a nonzero area). 
An example of such case for 02(0) = 0.5 is shown in Fig. HI where the solid curves represent 
the interface at various times t starting from t = 0.5 up to t = 2.0 with a time separation 
of At = 0.3 between successive curves, while the dashed curve represents the path traced 
by the tip j t = g^~ 1 (a2(t)). As one can see in this figure, for longer times the tip moves 
along a straight line in the z-plane, indicating that the initial asymmetry persists for all 
times. [In this and subsequent figures the trajectories of the tips and troughs are computed, 
for clarity, up to a final time that is slightly greater than the final time for the complete 
interfaces.] To generate the curves shown in Fig. H] we used a numerical scheme similar 



13 



-2 -10 12 

x 

FIG. 5: Symmetric growing interface with two tips. Here d,2 = d^ = — 1 and d$ = 0.5. The solid 
curves represent the interface at times from t = 0.5 to t = 1.7, with At = 0.3 between sucessive 
curves, whereas the dashed lines indicate the trajectories of the tips and the trough. 

to that described in Sec. [TT] for the case of multifingers, namely, we integrate the Loewner 
equation ( 12~4"|) backwards in time with terminal conditions g t = w, for w G (ai(t) , a^t)) , to 
get the corresponding points z = g on the interface. 

Next we consider the case in which the interface has two tips and one trough. In this 
case we have N = 5, and from ( 129]) and fl30|) one can easily write the parameters d\ and 
g?5 in terms of the growth factors c?2 and d± for the tips and d% for the trough. In Fig. [5] 
we show an example of a symmetrical growing interface (solid curves), where the dashed 
lines represent the trajectories of the tips and the trough. In this figure, we started with 
symmetrical initial conditions, namely, as(0) = — ai(0) = 1, 04(0) = — a2(0) = 0.5, and 
a 3(0) = 0, and chose the growth factors of the two tips to be the same, g?2 = = —1, so as 
to preserve the initial symmetry, with d% = 0.5. Note that the paths traced by the two tips 
and the trough are somewhat reminiscent of the trajectories observed in the case of three 



symmetrical fingers discussed in Sec. [TTJ see Fig. 2(a) On the other hand, if we start with 
a symmetric initial condition but one of the tips has a larger growth factor, then the fastest 
tip will move ahead of the other tip thus breaking the initial symmetry. For t — > 00 the 
slower tip and the trough will both tend to merge with the closest endpoint resulting in a 
growing interface with only one tip. An example of this case is shown in Fig. |5] where the 
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FIG. 6: Asymmetric growing interface with two tips. The initial conditions are the same as in 
Fig. [5]but here we have c?2 = —1, d% = 0.8, and = —0.5. 

tendency of the trough and the second tip to merge is clearly evident. The behavior seen 
in this figure mimics somewhat the competition between initial protuberances in the early 



12|. (In our model the competition is of course encoded 



stages of the fingering instability 
in the choice of the growth factors.) 

The generic Loewner equation ff24|) can also describe the problem of multiple growing 
interfaces. In this case each interface T l t , for i — 1, ...,n, where n is the number of distinct 
interfaces, will be mapped under gt(z) to a corresponding interval on the real axis in the 
w-plane. Similarly, each advanced interface T] +T is mapped by gt(z) to a polygonal curve 
in the w-plane. One can then readily convince oneself that the generic Loewner evolution 
defined by fl24"|) and ff28l) applies to this well, where N is now the total number of 

vertices corresponding to the sum of the number of vertices of each interface. 

In Fig. [7J we show numerical solutions for the case of two symmetrical interfaces with 
one tip each. Since each interface has three special points (two endpoints and one tip) we 
then have N = 6. Here the interfaces start to grow from the intervals [—3, —1] and [1,3], 
respectively, with the corresponding tips starting at symmetrical points a 4 (0) = —a 2 (0) = 2 
and having the same growth factors. The curves shown in Fig. [7J were calculated for d 2 = 
g?4 = —1, at times varying from t = 0.5 up t = 1.4 with a time interval At = 0.3 between 
successive curves. Notice that, as time goes by, the inner sides of the two interfaces move 
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FIG. 7: Loewner evolution for two symmetrical growing interfaces. The interfaces start respectively 
at [—3, —1] and [1, 3], and the growth factor is \d\ = 1 for each tip. The solid curves represent the 
interfaces at times varying from t = 0.5 up t = 1.4, with time intervals of At = 0.3, whereas the 
dashed line indicate the trajectories of the tips. 

towards one another leaving a narrow channel between them. For sufficiently large time, the 
width of such a channel becomes infinitesimally small so that for all practical purposes the 
resulting evolution will look like a single symmetrical growing interface. An evidence of this 
fact is clearly seen in Fig. [7] where the two tips (dashed lines) tend to "attract" each other 
and will eventually "merge." It is worth mentioning, however, that it is hard to integrate 
the Loewner equation past the latest time shown in Fig. [7J because the functions a^{t) and 
a$(t) become essentially identical (within the computer resolution), rendering the numerical 
integration very difficult after this point. 

In Fig. [8] we show for comparison the case in which the two tips have identical growth 
factors but where the interfaces now have different initial widths. We can see from this figure 
that the wider interface grows faster than the narrower one. This competition between the 
two growing interfaces resembles the so-called "shadowing effect" in fingering phenomena, 
whereby the longer fingers grow faster and hinder the growth of the shorter fingers in their 
vicinities 7]. [A similar effect is obtained if we start with interfaces of equal width but 
with one of the tips having a larger growth factor than the other one.] We have thus seen 
that in spite of its simplifying assumptions, chiefly among them the unbounded geometry 
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FIG. 8: Loewner evolution for two asymmetric growing surfaces. Here the growth factors and the 
initial width of first interface are the same as in Fig. [7J but the second interface starts to grow from 
the narrower interval [1.5, 2.5]. The different curves are for the same times as in Fig. [7J 

and the polygonal growth rule, our model exhibits certain dynamical features such as finger 
competition and screening that are qualitatively similar to what is commonly observed in 
Laplacian growth processes [ijj]. We will show next that it is possible, in principle, to 
extend the model to include more general growth rules that may describe more accurately 
the dynamics of specific growth processes. 

C. Loewner Domains 

Here we consider from a more formal viewpoint the Loewner evolution of a family of 
increasing hulls K t in the upper-half plane H. We denote by T t = K t H H the growing 
interface at time t and consider the mapping g t : M\K t — > H that maps T t to the interval 
[a(t),b(t)\ on the real axis. To specify the growth rule let us consider the infinitesimal hull 
K T t = g(K t+T \K t ) in the w-plane, where r is an infinitesimal time step. Note that in the 
growth models of Sec. IIII Al the hull K t;T corresponds to the interior of the polygon denned 
by g t (T t+T ) and the real axis; see Fig. [3j Here we consider more general growth problems 
in which K tjT is not necessarily bounded by a piecewise-linear curve in H. To be more 
precise, we consider the infinitesimal hull K r t = {w = x + iry, < y < h t (x)} made of 
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the set of points included between the real axis and the curve y = Th t (x), r < 1, with 
x real. We assume that h t (x) is a simple curve that is at least twice differentiable. Let 
us now denote by a*{t) < a(t) and 6*(t) > b{t) the endpoints of the curve h t (x), that is, 
h t (a*(t)) = h t (b*(t)) = 0, where \im T ^ a*(t) = a(t) and lim r ^ 6*(t) = b(t). Note that in 
the general case one has a* T {t) 7^ a(t) and b*(t) 7^ b(t), which allows for the endpoints Z\{t) 
and Z2(t) of the growing interface T t to move as time goes by. Thus, in contrast with the 
model discussed earlier, here the endpoints are not necessarily kept fixed. 

Let us now introduce a partition of the interval [a* (t), b*(t)]: a*(t) = ai(t) < a 2 (t) < 
■ ■ ■ < ajv(t) = b*(t), and let h t (x) the piecewise-linear approximation of h t (x) defined by 
this partition: h t (ai) = h(di), with h t (x), for x G [ai(t),a i+ i(t)}, being given by a straight 
line. It is clear that the evolution entailed by the polygonal curve h t (x) is described in terms 
of the Loewner equation ff24l) . If we now take the limit N — > 00, we then get the following 
evolution equation: 

g t (z) = / K t (x)[g t (z) - x] hx[g t {z) - x}dx, (32) 

Ja(t) 

where 

Kt(x) = h"(x), (33) 

with prime denoting derivative with respect to the argument x. Since the shape of the 
growing domain is fully encoded in the map gt(z), equation ( |32l) specifies the growth model 
once the density Kt(x) is known. Note that the class of interface growth models described 
in Sec. IIII Al corresponds to the case when the density is a finite sum of Dirac 5 peaks. 
More generally, if we introduce a signed measure /it{x) satisfying the conditions 



dfit(x) = 0, (34) 

/ xd[i t (x) = 0, (35) 
Jr 

which are the analog of ( 1291) and ( |30l) . we can then define the Loewner evolution generated 
by this measure by the following equation 

9t(z) = / \gt(z) - x] \a[g t {z) - x]dfi t (x). (36) 
Jr 

Here the time evolution of the measure m(x) encodes the growth rule. We refer to ( 136]) 
as "Loewner domains" in contrast to the "Loewner chains" discussed in Sec. [Til Although 
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Loewner domains can be regarded as equivalent to Loewner chains — indeed, the densities 
Pt(x) and Kt(x) are obviously related to one another — , the former description seems to be 
more convenient for the study of growing interfaces (which encircle domains with nonzero 
area). From a more practical viewpoint, however, the difficult task that remains in the 
case of Loewner domains is to find the appropriate measure for physically relevant growth 
problems. In this context, an interesting open question concerns the possibility of generating 
rough (i.e., fractal) interfaces by means of a Loewner evolution driven by random measures, 
in analogy with the stochastic Loewner equation where random curves are produced % \5\. 



IV. SUMMARY AND CONCLUSIONS 



We have considered the problem of two-dimensional Laplacian growth in the context of the 
Loewner evolutions. We started by revisiting the problem of fingered growth and reported 
a novel exact solution of the chordal Loewner equation for the case of a symmetrical three- 
finger configuration. We then discussed a more general class of growth models in which 
an interface grows from the real axis into the upper-half plane. Assuming that the growth 
rule is specified in terms of a polygonal domain in the complex-potential plane and making 
appropriate use of the Schwarz-Christoffel transformation, we derived the corresponding 
Loewner equation for the problem. Several examples were explicitly discussed, such as the 
case of a single "extended finger," an interface with two tips, and the case of two competing 
interfaces. Although our model does not allow for a direct comparison with experiments, we 
have argued that its dynamics is reminiscent of certain typical behaviors seen in Laplacian 
growth. We have also briefly discussed a possible extension of our Loewner equation to 
include the case in which the growth rule is specified in terms of a generic (not necessarily 
polygonal) curve. More generally, we have introduced the notion of "Loewner domains" in 
which the growth rule is described in terms of a time evolving signed measure. This opens 
up the possibility for studying other interesting (and difficult) problems, such as the growth 
of fractal interfaces, within the context of stochastic Loewner evolutions. The extension to 
the radial geometry of the class of growth models discussed here is another interesting open 
problem. 
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